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ABSTRACT 

The differing M^-L relations presented in McLure & Dunlop, Marconi & Hunt 
and Erwin ct al. have been investigated. A number of issues have been identified and 
addressed in each of these studies, including but not limited to: the removal of a de- 
pendency on the Hubble constant; a correction for dust attenuation in the bulges of 
disc galaxies; the identification of lenticular galaxies previously treated as elliptical 
galaxies; and application of the same (^I^Q regression analysis. These adjustments 
result in relations which now predict similar black hole masses. The optimal if-band 
relation is log(Af b h/Af ) = -0.37(±0.04) [M K + 24] + 8.29(±0.08), with a total (not 
intrinsic) scatter in logAfbh equal to 0.33 dex. This level of scatter is similar to the 
value of 0.34 dex from the Afbh-c relation of Trcmainc et al. and compares favourably 
with the value of 0.31 dex from the M^-n relation of Graham & Driver. Using dif- 
ferent photometric data, consistent relations in the B- and i?-band arc also provided, 
although we do note that the small (N = 13) i?-band sample used by Erwin et al. is 
found here to have a slope of —0.30 ± 0.06 and a total scatter of 0.31 dex. Perform- 
ing a symmetrical regression on the larger if-band sample gives a slope of ~ —0.40, 
implying Mbh oc L 100 . Implications for galaxy-black hole coevolution, in terms of dry 
mergers, are briefly discussed, as are the predictions for intermediate mass black holes. 
Finally, as previously noted by Tundo et al., a potential bias in the galaxy sample used 
to define the M^-L relations is shown and a corrective formula provided. 
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1 INTRODUCTION 

The mass of a spheroid's supermassive black hole (SMBH), 
denoted by Mbh, is known to correlate with several physical 
properties of the spheroid, by which we mean either an el- 
liptical galaxy or the bulge of a disc galaxy. At first it was 
thought that the velocity dispersion, a, of the bulge was 
the primary driving mechanism (Ferrarese & Merritt 2000; 
Gebhardt et al. 2000). However, in the following year the 
radial concentration of stars in the spheroids was observed 
to correlate just as well with the SMBH mass and yielded 
a relation with the same small degree of scatter (Graham 
et al. 2001). While the luminosity of the stars initially ap- 
peared to provide a weaker relation, it is now known that this 
was predominantly due to rough estimates of the spheroid's 
luminosity. Performing a Sersic (1963) -R 1//n -bulge plus ex- 
ponential disc decomposition of the galaxy light, Marconi & 
Hunt (2003) and Erwin et al. (2004) revealed that the total 
scatter about the near-infrared and optical Mbh-L relations 
was similarly small at ~0.35 dex. Most recently, Graham 
& Driver (2007a, their Section 6) have predicted that the 
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central stellar density, prior to core-depletion of the host 
spheroid, may also be intimately connected with the SMBH 
mass. 

These relationships are important for two main reasons. 
First, they provide an easy means to predict SMBH masses 
in thousands of galaxies for which direct measurements of 
the SMBH mass is not possible (e.g., Yu & Tremaine 2002; 
Marconi et al. 2004; Shankar et al. 2004; Graham et al. 2007 
and references therein). The other reason is that they are a 
clue to the driving physical processes at work in galaxies. 

However, as with the controversy over the slope of the 
Mbh-f relation (Merritt & Ferrarese 2001a; Tremaine et al. 
2002; Novak et al. 2006), the slope of the M hh -L relation 
is not yet agreed upon. While McLure & Dunlop report a 
value of -0.50 ± 0.05, the data in Erwin et al. (2004) has 
a slope of —0.25 ± 0.05. Moreover, the differing relations 
from different studies tend not to predict the same SMBH 
mass. For example, when Lk = 10 10 L K q (10 11 L k q), the 
if-band expression in McLure & Dunlop (2004, their equa- 
tion 1) predicts SMBH masses which are three (two) times 
less than those predicted by the if-band Mbh-L relation in 
Marconi & Hunt (2003). 

In this paper we re-investigate the past Mbh-L rela- 
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tions and address a number of ways in which they can be 
updated. In Section 2 we present an early data set used to 
construct an Mbb-L relation. We use this data to illustrate 
how we define the regression analysis that we shall adopt in 
this paper. Section 3 starts with studies which avoided the 
bulge/disc separation issue by excluding disc galaxies, while 
Section 4 explores those studies which used elliptical, lentic- 
ular and spiral galaxies. We obtain new B-, R- and A'-band 
relations which are consistent with each other and suitable 
for predicting black hole masses in other galaxies. In Sec- 
tion 5 we discuss how intermediate mass black holes factor 
in, and briefly mention some of the implications of these new 
relations for the coevolution of spheroids and SMBHs, 

Finally, in an appendix we provide a re-derivation of 
the Mbb-L relation taking into allowance a potential bias in 
the galaxy sample. 



which has a total rms scatter of 0.56 dex in the log Mbh 
direction and an intrinsic scatter (assuming a magnitude 
error of 0.3 mag) of 0.46t^| dex. 

This level of scatter is unpleasantly high and resulted 
in the Mbb-L relation taking second place to the Mbh- 
cr (and Mbh-w) relation. However, it has since been re- 
alised/shown that this level of scatter was a consequence 
of a poor bulge/disc separation (Marconi & Hunt 2003; Er- 
win et al. 2004) and the use of systems in which the SMBH's 
sphere of influence was not well resolved (e.g. Merritt & Fer- 
rarese 2001c; Ferrarese & Ford 2005). The following Section 
avoids the issue of the bulge/disc separation by dealing with 
a study that used elliptical galaxies. Section [3] effectively 
tackles this issue by using studies in which a Sersic-bulgq 2 ] 
plus exponential-disc decomposition of the galaxy's stellar 
light has been performed. 



2 CONSTRUCTION OF THE Mbk-L 
RELATION 

2.1 Kormendy & Gebhardt (2001) 

As noted by Kormendy & Gebhardt (2001, see their sec- 
tion 5) SMBHs seem to be associated with the dynami- 
cally hot, spheroidal component of a galaxy. Those authors 
therefore presented a B-band relation based on the bulge, 
rather than total, magnitudes from Faber et al. (1997). Us- 
ing M B q = 5.47 mag, their Mbh-L relation can be written 
as 



log(M bh /M ) = -0.43[Af s + 19.5] + 7.1 



(1) 



Given our objective is to construct the optimal rela- 
tion (y = a + bx) for estimating the SMBH mass (y) from 
the magnitude (x) of a galaxy, we adopt the method of re- 
gression analysis given in Tremaine et al. (20020 We allow 
for intrinsic scatter (in the y-direction) , which we denote 
by the term e, and also for measurement errors on the N 
pairs of observables Xi and yi, which we denote 5xi and Syi. 
Tremaine et al.'s (2002) modified version of the routine FI- 
TEXY (Press et al. 1992, their Section 15.3) minimises the 
quantity 



X 



N 

£ 



bxA 



5yS + b 2 5x^ +e 2 



(2) 



The intrinsic scatter e is solved for by repeating the fit until 
X 2 /(N — 2) equals 1. The uncertainty on e is obtained when 
the reduced chi-squared value, x 2 /{N — 2), equals lzb-^/2/A' . 
To achieve a minimisation in the ^-direction, one simply re- 
places the e 2 term in the denominator of equation [2] with 
6 2 e 2 . A symmetrical regression is therefore some kind of av- 
erage of these two regressions (Novak et al. 2006) . 

Application of equation [2] to the data from Kormendy 
& Gebhardt yields the relation 

log(M bh /M ) = -0.38(±0.06)[M s +19.5]+8.00(±0.09), (3) 



3 ELLIPTICAL GALAXIES 
3.1 McLure & Dunlop (2002) 

Given that the mass of a SMBH is known to correlate with 
the properties of the host spheroid, rather than the host 
galaxy, McLure & Dunlop (2002) presented the Mbh-L re- 
lation after excluding (the bulk of the) disc galaxie^f]. Their 
expression (for a predominantly elliptical galaxy sample) is 



log(M bh /M ) = -0.50(±0.05)M fl - 2.91(±1.23). 



(4) 



Here we re-derive this expression after implementing the fol- 
lowing alterations. 

• Measurements of black hole mass depend linearly on the 
distance to each galaxy. Distances for 16 of the 18 galax- 
ies used by McLure & Dunlop have been obtained using 
surface brightness fluctuations (Tonry et al. 2001) and so 
their SMBH masses are independent of the Hubble constant. 
However, in converting the apparent magnitude of these 18 
galaxies into absolute magnitudes, McLure & Dunlop used 
the galaxy redshift and a Hubble constant Ho = 50 km s _1 
Mpc - . This was done because of their comparison with a 
sample of AGN for which Ho = 50 km s~ Mpc -1 had been 
used. A more consistent approach would involve the use of 
the //o-independent distances for the non-AGN galaxies to 
determine their absolute magnitudes, and a Hubble constant 
of 73 km s" 1 Mpc" 1 (Blakeslee et al. 2002) for the AGN 
sample. Here we use these //o-independent distances for both 
the SMBH masses and the spheroid magnitudes of our (local 
non-AGN) galaxies. 

• Nine of the 18 7?-band magnitudes which McLure & 
Dunlop used were (B — R c = 1.57)-adjusted B-band mag- 
nitudes from Faber et al. (1997) — which themselves came 
from the "Seven Samurai" data (Faber et al. 1989) and/or 
the third reference galaxy catalogue (RC3, de Vaucouleurs et 
al. 1991). The other nine i?-band magnitudes are reported to 
be (V — R — 0.61)-adjusted F-band magnitudes from Mer- 
ritt & Ferrarese (2001b). Although this paper has no F-band 



1 This is not a symmetrical method of regression, but rather one 
which minimises the scatter in the y-direction (see Novak et al. 
2006 and Graham & Driver 2007a). The criticism in Tundo et al. 
(2007, their Section 2.2) of the Tremaine et al. (2002) method is 
therefore misplaced. 



2 A modern review of the Sersic model can be found in Graham 
& Driver (2005). 

3 Bettoni et al. (2003) performed the same task, obtaining a con- 
sistent relation. 
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magnitudes, Merritt & Ferrarese (2001a) has V-band mag- 
nitudes from Faber et al. (1989) for 4 of these 9 galaxies. We 
therefore use the Faber et al. (1989) B-band apparent mag- 
nitudes (and the related absolute magnitudes in Tremaine 
et al. 2002). The Faber et al. (1989) magnitudes are derived 
from photoelectric aperture growth curves. While the mag- 
nitude quality indicators from that paper suggest an error 
of < 0.15 or 0.30 mag, depending on the galaxy, we have 
adopted the upper value for all galaxies in our regression 
analysis. 

• We have removed NGC 4564 and NGC 2778 which are 
not elliptical galaxies but SO galaxie^f] (see Graham & Driver 
2007b) whose total galaxy luminosities would have biased 
the previous relation. We have also excluded the peculiar 
elliptical galaxy IC 1459 due to uncertainty on its SMBH 
mass. While the stellar dynamics of its core suggest a SMBH 
mass of 2.6x 10 9 Mq, the gas dynamics reveal the mass could 
be as low as 3.5 x 10 8 Af Q (Cappellari et al. 2002). We have 
treated NGC 221 as an SO galaxy according to the B/T flux 
ratio in Graham (2002). 

• We have included NGC 1399 and NGC 5845 for which 
accurate SMBH masses have since become available, giving 
a total of 17 galaxies. 

• Due to our desire to obtain a relation for predicting ac- 
curate SMBH masses using the magnitudes of other galax- 
ies, we perform a non-symmetrical regression which results 
in the smallest degree of scatter in the logMbh direction. 
That is, we apply equation [5] 

Table [1] provides distances, updated SMBH masses 
(based on these distances) and absolute magnitudes (again 
based on these distances) for this set of galaxies. Using the 
distances and apparent bulge magnitudes from Table [TJ and 
applying equation [5] we obtain the relation 

log(Af b h/M ) = -0.42(±0.06)[M S + 20] + 8.32(±0.10), (5) 

with an absolute scatter in log Mbh of 0.36 dex. This relation 
can be seen in Figure[T] Using the absolute magnitudes given 
in Table[TJ rather than the apparent magnitudes, one obtains 

log(Af b h/M ) = (-0.36±0.06)[M B +20] + (8.33±0.10).(6) 

These equations have an order of magnitude less uncer- 
tainty on the intercept term than the uncertainty given in 
equation 2] They also have a shallower slope. 

3.1.1 Independence of Hq 

As can be seen in Table [TJ two of the galaxies used to con- 
struct the above Mhh-L relation have magnitudes and black 
hole masses that depend on the Hubble constant. We have 
explored whether or not equation [S] and [5] would be signifi- 
cantly different if (i) these two galaxies were excluded and 
(ii) if we had used a Hubble constant of 50 or 100, rather 
than 73 km s _1 Mpc -1 for these two galaxies (NGC 6251 
and NGC 7052). In all cases the slope and intercept var- 
ied by no more than 0.02 and 0.03, respectively. What this 
means is that the M^-L relations given above, and in the 

4 Lenticular galaxies have typical bulgc-to-total luminosity ratios 
of 1/4 (e.g., Balcells, Graham & Peletier 2004; Laurikainen, Salo 
& Buta 2005). 



Table 1. Revised (see Section 13.11 1 sample of elliptical galaxies 
from McLure & Dunlop (2002). Distances are taken from Tonry 
et al. (2001, their table 1), except for NGC 6251 (i' C mb=7382 km 
B -1 , Wcgncr ct al. 2003) and NGC 7052 (i> CMB =4411 km s -1 , 
Wegner et al. 2003). These two galaxies are not listed in Tonry 
et al. (2001) and a Hubble constant of Hq = 73 km s — 1 Mpc -1 
(Blakeslee et al. 2002; Spergel ct al. 2006) has been used. Unless 
noted otherwise, the apparent B-band magnitudes, mj, are the 
_By magnitudes from Faber et al. (1989, 1997). Unless noted oth- 
erwise, the absolute B-band magnitudes, Mg, have come from 
Tremaine et al. (2002) with modifications for NGC 6251 and 
NGC 7052 such that we adjusted the Tremaine et al. values 
(-21.81 and -21.31 mag) to our adopted distance. The SMBH 
masses are also from the compilation in Tremaine ct al. (2002), 
except for NGC 821 (Richstone et al. 2007), NGC 3379 (Gebhardt 
et al. 2000; see also Shapiro et al. 2006), and NGC 4486 (Mac- 
chetto et al. 1997). Our sample includes two additional galaxies 
not used in Tremaine et al. The SMBH mass for NGC 1399 is 
from Houghton et al. (2006) and the mass for NGC 4374 is from 
Macicjewski & Binney (2001, with updated errors taken from 
Kormendy & Gebhardt 2001). 
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a Reduced using a B/T ratio of 0.62 (Graham 2002). 

" R-band bulge magnitude from Graham (2002) with a B — R = 

1.84 mag adjustment (Lugger et al. 1992. 

c SMBH sphere of influence not resolved. 

^ Derived from the apparent magnitude. 

e SMBH sphere of influence not resolved. 

/ Derived from the apparent magnitude. 

9 This SMBH mass is based on M.E. Kaiser (2001, in prep.). 
h RC3 B(m B ) value. (No value in Faber et al. 1989.) 
• RC3 B(m° B ) value. (No value in Faber ct al. 1989.) 



rest of this Section, are effectively independent of the Hubble 
constant. 



3.1.2 An R-band Mbh-L relation 

As already noted, the i?-band relation from McLure & Dun- 
lop (2002) was based on the Faber et al. (1989) B-band mag- 
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Figure 1. Correlation between an elliptical galaxy's supcrmassive black hole mass and the apparent S-band magnitudes listed in 
Table [T] The regression line shown in the left panel was obtained using equation [5] NGC 221 (M32) is plotted with a circle around it. 
The middle panel shows the A\ 2 = 1.0 and 2.3 boundaries around the optimal intercept, a = 8.32, and slope, b = —0.42. The projection 
of the Ax 2 = 1.0 ellipse onto the vertical and horizontal axis gives the 1-a uncertainties Sa and 8b, respectively. The Ax 2 = 2.3 ellipse 
denotes the 1-c two-dimensional confidence region. This latter quantity has been mapped into the right panel, and is traced by the two 
solid curves. The dashed lines in this panel are the (more commonly used) approximations obtained using (a±5o) and (6 ±56). The two 
confidence regions agree well, although the region traced by the dashed lines is, as expected, smaller. 



nitudes and an B - R c colouiQ equal to 1.57 mag. Adopting 
this same colour term, except for NGC 221 which has a B — R 
colour of 1.84 (Lugger et al. 1992), and using the absolute 
magnitudes given in Table [T] one obtains the relation 

log(M b h/M ) = -0.38(±0.06)[M fl + 21] + 8.11(±0.11). (7) 

For a dynamically hot spheroid with Mr — —21 mag, 
McLure & Dunlop's relation (equation [4| gives values of 
log Mbh which are 0.52 dex less massive, i.e. roughly a factor 
of three less massive, and considerably offset from our refined 
zero-point at Mr = —21 mag, given by 8.11±0.11. However 
it should be noted that most of this offset (at Mr = —21 
mag) is simply due to the prior use of Ho = 50 km s _1 
Mpc -1 when deriving the absolute magnitudes used in equa- 
tion|3] If a value of 73 km s _1 Mpc -1 is used there, then the 
masses are now only 0.11 dex smaller, which is consistent 
with the 0.11 dex uncertainty on the new zero-point. 



3.1.3 A K-band M hh -L relation 

We also update the i^-band relation presented in McLure 
& Dunlop (2004, their equation 1), which was derived there 
using an R — K colour of 2.7 mag applied to their J?-band 
relation from McLure & Dunlop (2002). When McLure & 
Dunlop (2004) converted McLure & Dunlop's (2002) M bh -L 
relation, they also applied a Hubble conversion of log(50/70) 
to the black hole masses and 5 log(50/70) to the magnitudes. 
This is not applied here because equation [7] is, for prac- 
tical purposes, independent of the Hubble constant (Sec- 
tion [3TLT}. Adopting an R — K colour of 2.6 mag (Buzzoni 
2005), equation [7] becomes 

log(Af b h/M ) = -0.38(±0.06)[A/ A - +24] + 8.26(±0.11).(8) 

Using here an absolute JT-band magnitude for the Sun of 
M Xj q = 3.28, this can alternatively be expressed as 

5 From here on the subscript c shall be dropped from the term 
Rc- 



log(M bh /M Q ) = 0.95(±0.15) log .J^f +8.26(±0.11), (9) 

where Lk,s P Ii/L k q is the K-band luminosity of the 
spheroid component of the galaxy (i.e., the bulge or the el- 
liptical galaxy itself) in solar units. 



4 INCLUDING THE DISC GALAXIES 

More recent studies have included both elliptical galaxies 
and the bulges of disc galaxies. In this section we explore 
the Mbh-i expressions obtained with both sets of objects. 

4.1 Erwin, Graham & Caon (2002) 

In 2002 Erwin, Graham & Caon presented a relation be- 
tween black hole mass and host spheroid magnitude for a 
sample of 13 galaxies (8 elliptical and 5 disc galaxies). The 
study was presented at the Carnegie Observatories Astro- 
physics conference "Coevolution of Black Holes and Galax- 
ies" and posted to astro-ph that same year. For the disc 
galaxies, the bulge magnitudes were obtained from an R}^ n - 
bulge plus exponential-disc decomposition (Graham et al. 
2001). This was the first study, albeit small in number, to 
show that the inclusion of bulge galaxies resulted in a total 
scatter of 0.35 dex in log Mbh, considerably less than past 
reports of ~0.6 dex and comparable to the scatter in the 
Mbh-o" and Mhh-n relations. 

Using the data points shown in Figure 1.4 of Erwin et al. 
(2004), along with a 0.3 mag uncertainty on the magnitudes, 
we have performed a regression of log Mbh against these it- 
band magnitudes (i.e., applied equation [2} . Doing so gives 

log(M b h/M ) = -0.24 ± QM[Mr + 21] + 8.01 ± 0.11, (10) 

with a total scatter of 0.35 dex. This is shown in Figure [2ji. 

This relation has a noticeably shallower slope than the 
value in equation [7] of — 0.38(±0.06) obtained using our 
refined elliptical galaxy sample from McLure & Dunlop 
(2002). When Mr = —24 mag, equation 1101 gives masses 
of ~ 5 x 10 s Mq — three times smaller than that obtained 
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Figure 2. The long solid lino shows the M^-L relation (equa- 
tion 1101 and 1121 using the 8 elliptical galaxies (filled circles) plus 
the bulges of 5 lenticular /spiral galaxies (open circles) given in 
Erwin et al. (2004). The short lines emanating from the elliptical 
galaxy data points in panel a) show the location of these galaxies 
as used in Section 13.11 to obtain equation [7] (shown by the long 
dashed line). Panel a) shows the data points as seen in Erwin et al. 
(2004; their Figure 1.3). Panel b) shows the data points as given 
in Table [2] which uses slightly updated SMBH masses and, most 
importantly, dust-corrected bulge magnitudes (see Section 1 1.1. 11 1. 
The dotted line has a slope of —0.4. 



using equation [7] When Mr = —27 mag, equation [TOl gives 
masses of ~ 3 x 10 9 AfQ. In this regard, equation \W\ does not 
conflict with the Mhh-c relation. That is, compared to the 
Mbh-o" relation, equation [10] does not predict significantly 
larger SMBH masses at the high mass end. However, even 
if the Sersic bulge (and elliptical galaxy) magnitudes used 
in Erwin et al. are reliable, the low number of points (only 
13) may make equation \W\ prone to statistical fluctuations 
in the selected sample. 

In Figure[2]we show the data points used by Erwin et al. 
(2004). Also shown, for the elliptical galaxies, is their loca- 
tion using the values adopted in the previous Section where 
a slope of —0.38 was obtained. The only data point which 
has shifted significantly is NGC 821. The reason is because 
the previous section used the SMBH mass from Richstone 
et al. (2007, 8.5 x 10 7 A/ Q ), while Erwin et al. used the (at 
the time unpublished) value of 3.7 x 10 7 Mq (Tremaine et 
al. 2002). 



4.1.1 Dust 

While, in general, elliptical galaxies and the bulges of disc 
galaxies are considered not to have dust, this does not mean 
that the dust in the discs of disc galaxies does not influence 
the emergent flux from the stars in the bulge. To illustrate 
this, consider an infinitely thin and optically opaque dust 
sheet running through the disc of a galaxy. Obviously, near 
edge-on orientations aside, one will only see the bulge (and 
disc) stars on the near-side of the disc. That is, one will only 
see half of the optical flux from the bulge, it will therefore 
be observed to be 0.75 mag fainter than it actually is. In 
reality, dust discs have a certain thickness, relative to the 
vertical scale-height of the stars in the disc. This results in 
a dust correction that depends on the inclination of the disc 
relative to our line of sight. 

Here we modify the B-band inclination-attenuation cor- 
rection for bulge magnitudes given in Driver et al. (2007, 
their equation 3). Specifically, we use the expression 



in which the 0.84 term provides the face-on attenuation- 
correction to the observed magnitude of the bulge, M b s , 
while the latter part of the expression provides the 
inclination-dependent component of the correction. Equa- 
tion [11] only differs from the B-band relation due to our 
(simplistic) use of a 0.6 multiplier. This multiplicative factor 
stems from our knowledge of the disc extinction/attenuation 
in the B- and .R-bands (e.g. Tully & Verheijen 1997; Tully 
et al. 1998). 

The five bulge magnitudes listed in Tableware corrected 
for dust using equation llll The values plotted in Erwin et al. 
(2004), and the left panel of Figure [2] are some 0.55 to 0.80 
mag fainter. Re-performing the regression analysis, with the 
dust- corrected magnitudes, and using the Richstone et al. 
(2007) SMBH mass for NGC 821, we obtain 

log(M bh /M ) = -0.30(±0.06)[Mr+21] + 7.96(±0.10), (12) 

with a scatter of 0.31 dex in logAfbh- This is shown in Fig- 
ure [2}). 

We have repeated this regression after the jackknife re- 
moval of NGC 2787, a galaxy with an inner disc twice as 
luminous as its bulge (Erwin et al. 2003), and the faintest 
object in our sample. Doing so results in a slope and inter- 
cept of — 0.35T0.05 and 7.90T0.09, respectively, and a total 
scatter of 0.29 dex. 

To better gauge the uncertainty on the slope in equa- 
tion 1121 we have also used a bootstrap sampling of 13 
data points (i.e. sampling with replacement from the origi- 
nal sample) and 1000 such Monte Carlo samples. The me- 
dian ±2<7 of the resultant distribution of 1000 slopes was 
-0.31 ± 0.11. For the intercept we obtained 7.961^. 

4.2 Marconi & Hunt (2003) 

4-2.1 A K-band M^-L relation 

Using images from 2MAS^ Marconi & Hunt (2003) ob- 
tained K-hanA magnitudes for 27 galaxies which had di- 
rect and reliable SMBH mass measurements. Recognising 
the need for bulge magnitudes rather than galaxy magni- 
tudes, they also performed an _R 1 ^ n -bulge plus exponential- 
disc decomposition of the disc galaxies in their sample. With 
a sample size twice that used in Erwin et al. (2004), they con- 
firmed Erwin et al.'s claim for a scatter of around 0.35 dex 
in the Mbh-L relation for elliptical galaxies and the bulges 
of disc galaxies. In what follows we describe the manner and 
reason for why we have tweaked the masses and magnitudes 
for some of their galaxy sample. 

• Due to the shallow nature of the 2MASS images (only 8 
seconds) the presence of a disc in NGC 27783 and NGC 4564 
were missed. NGC 221 (M32) was also treated as a pure 
elliptical galaxy. From the i?-band bulge/disc decomposition 
of these three galaxies in Graham & Driver (2007a) and 



6 Two micron all sky survey: Jarrett et al. (2000). 

7 Because we used Marconi & Hunt's preferred 27 "Group 1" 
galaxies (see their Table 1), NGC 2778 does not actually factor 
into the analysis here because it is tabulated as a "Group 2" 
galaxy. 
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Table 2. Revised sample of elliptical and disc galaxies from 
Erwin et al. (2004). A galaxy Type 'E' denotes an elliptical galaxy 
while a Type 'S' denotes either an SO galaxy or a spiral galaxy. 
Distances are taken from Tonry et al. (2001, their table 1), except 
for NGC 7052 (« CMB =4411 km s -1 , Wegner et al. 2003) which is 
not listed in Tonry et al. A Hubble constant of Ho = 73 km s _1 
Mpc — 1 has been used for this galaxy. The SMBH masses are from 
the compilation in Tremaine et al. (2002). except for NGC 821 
(Richstone et al. 2007) and NGC 3379 (Gebhardt et al. 2000). Our 
sample includes one additional galaxy not used in Tremaine et al. 
The SMBH mass for NGC 4374 is from Maciejewski & Binney 
(2001, with updated errors taken from Kormendy & Gebhardt 
2001). Disc inclinations, i.e. b/a axis ratios of the outer isophotcs, 
have come from: NGC 2778 (Rix et al. 1999); NGC 2787 (Erwin 
et al. 2003); NGC 3384 (Faber et al. 1997); NGC 4564 (Faber 
et al. 1997, see also Graham & Driver 2007a); and NGC 7457 
(Chapelon et al. 1999). The magnitudes are from Erwin et al. 
(2004), but corrected for dust attenuation using equation 1111 



Galaxy 



NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 



0821 
3377 
3379 
4261 
4374 
4473 
5845 
7052 
2778 
2787 
3384 
4564 
7457 



Type 

E 
E 
E 
E 
E 
E 
E 
E 
S 
S 

s 
s 
s 



Dist. 

[Mpc] 
24.1 
11.2 
10.6 
31.6 
18.4 
15.7 
25.9 

22.9 
7.5 
11.6 
15.0 
13.2 



b/a M R 



0.72 
0.57 
0.45 
0.45 
0.59 



[mag] 
-22.10 
-21.27 
-21.54 
-23.33 
-23.00 
-20.82 
-20.55 
-23.57 
-18.74 
-18.25 
-18.95 
-19.88 
-18.60 



[1O 8 M ] 

-+0^35 
-0.35 
-1+0.9 
-0.1 
;+0.73 
-0.73 
rl.O 
-1.1 
,+3.46 
-1.83 
H).40 
-0.79 
-,+0.4 
-1.4 
-,+2.4 
-1.3 
,+0.08 
-0.09 
1+0.04 
-0.05 
,+0.01 
-0.02 
hO.03 
-0.08 
-+0.011 
-0.014 



0.85] 
l.OOj 
1.35j 
5.2o] 
4.64^ 

l.ioj 

2.40j 
3.40^ 
0.14^ 
0.41^ 
0.16j 
0.56^ 
0.035^ 



Graham (2002), the bulge-to-total flux ratios are 0.21, 0.24 
and 0.62 respectively. We have used these ratios to obtain 
the bulge magnitudes for these galaxies. In passing we note 
that lenticular galaxies typically have B/T ratios of ~ 0.25± 
0.10 (Balcells et al. 2004; Laurikainen et al. 2005), the higher 
B/T value for M32 is to be expected if this is a partially disc- 
stripped lenticular galaxy (see Bekki et al. 2001; Graham 
2002). 

• We have excluded IC 1459 due to the order of magni- 
tude uncertainty on its black hole mass (Cappellari et al. 
2002), reducing our sample size from 27 to 26 galaxies (see 
Table EJ. 

• We have updated the SMBH mass and its associated 
uncertainty for NGC 5252 using the now published result in 
Capetti et al. (2005) together with a distance of 94.4 Mpc 
(slightly different to the value of 96.8 Mpc in Marconi & 
Hunt, and obtained using a recession velocity of 6888 km 
s _1 and H = 73 km s _1 Mpc -1 ) to give M hh = 0.97tJ;t§ x 
10 9 Mq. While this distance and mass is only 2.5 per cent 
smaller than that given in Marconi & Hunt, the uncertainty 
on the mass is 2-3 times larger. 

• We have also slightly modified the SMBH mass for 
Cygnus A by using the value M bh = (2.6 ± 0.7) x 1O 9 M . 
This was obtained from the mass in Tadhunter et al. (2003) 
after using a redshift of 0.056 and adopting Ho = 73 km s _1 
Mpc - together with fl m = 0.3 and FIa = 0.7. This gave a 
luminosity distance of 240 Mpc and an angular distance of 
215 Mpc (c.f. 207 Mpc in Tadhunter et al.). This mass is 



still consistent with the value of (2.9 ± 0.7) x 10 9 M@ used 
in Marconi & Hunt (2003). 

• In addition to NGC 5252 and Cygnus A, Marconi & 
Hunt's sample includes a further four galaxies not listed in 
Tremaine et al. (2002). The SMBH mass we have used for 
NGC 3031 (M81) is from Devereux et al. (2003), the mass for 
NGC 5128 (Cen A) has come from Marconi et al. (2001) and 
the mass for NGC 4594 is from Kormendy et al. (1988). The 
mass for NGC 4374 (M84) has been taken from Maciejewski 
& Binney (2001, with updated errors taken from Kormendy 
& Gebhardt 2001). These masses are the same as used by 
Marconi & Hunt, except for NGC 4374. 

• For the remaining 20 galaxies we have adopted the 
SMBH masses given in Tremaine et al. (2002), with only 
the following three exceptions: NGC 3379 (Gebhardt et al. 
2000; see also Shapiro et al. 2006), NGC 4486 (M87, Mac- 
chetto et al. 1997) and NGC 3115 (Emsellem, Dejonghe, & 
Bacon 1999). Our mass for NGC 3379 is 35 per cent larger 
than that used by Marconi & Hunt, while the difference in 
mass for the other two galaxies is only ~1 per cent from 
that used by Marconi & Hunt. 

• We also note that Marconi & Hunt used a distance of 
107 Mpc for NGC 6251, slightly greater than the value of 101 
Mpc which we adopted in the previous section and which we 
use here for consistency. This results in our reduction of the 
SMBH mass for this galaxy by the fraction 101/107, and a 
dimming of the absolute magnitude by 5 log(107/101). 

Table [3] presents our updated and modified data set 
from Marconi & Hunt's Group 1 galaxies. Finally, we again 
note that we perform a non-symmetrical regression analysis, 
as is given in equation [5] Marconi & Hunt used the symmet- 
rical bisector linear regression algorithm of Akritas & Ber- 
shady (1996) when reporting their optimal relation. As can 
be seen in their Figure 1, this results in a slightly steeper 
slope than obtained with an ordinary (non-symmetrical) 
least-squares fit to the data (see also Section \5. 2 1 of this pa- 
per). Marconi (2007, priv. comm.) reports that the slope of 
their data using equation [2] is —0.41 ± 0.04. 

Figure [3] shows the 7f-band Mhh-L relation derived us- 
ing our slightly updated Marconi & Hunt data set. One can 
also see, via the short lines, how each data point has moved 
from its previous location as given by Marconi & Hunt. Per- 
forming a linear regression (using equation [5]) , we obtain 

log(A/ bh /Af ) = -0.39(±0.05)[Mjf +24]+8.24(±0.08), (13) 

with a total scatter of 0.35 dex in log Mbh- This equation 
has a slightly shallower slope than the expression given by 
Marconi & Hunt (2003, their Table 2): log(M bh /M ) = 
1.13(±O.12)[log(LK,buigc/iK, ) - 10.9] + 8.21(±0.07), or 
equivalent lj|f|, using M K q = 3.28 mag, log(Af bh /M ) = 
-0.45(±0.05)[M J f,bui gc + 23.97] + 8.21(±0.07). However, as 
noted above, applying the same regression analysis to the 
data in Marconi & Hunt yields a consistent result. 



8 Dong & De Robertis (2006) obtain the same slope (—0.45) upon 
excluding the disc galaxies, but with a large uncertainty on the 
intercept. 
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Table 3. Updated sample of elliptical and disc galaxies from 
Marconi & Hunt's (2003) Group 1 galaxies. When available, dis- 
tances from Tonry et al. (2001) have been used. The excep- 
tions are NGC 5252, NGC 6251 and Cygnus A, as noted in Sec- 
tion [472TT] The magnitudes have come from Table 1 in Marconi & 
Hunt, adjusted here if a different distance was adopted (3 galax- 
ies) or if we assigned a bulge-to-disc ratio not used by Marconi 
& Hunt (NGC 221 and NGC 4564). The source of the black hole 
masses is also provided in Section l4.2.1l 



Galaxy 

NGC 0221 
NGC 1023 
NGC 2787 
NGC 3031 
NGC 3115 
NGC 3245 
NGC 3377 
NGC 3379 
NGC 3384 
NGC 3608 
NGC 4258 
NGC 4261 
NGC 4291 
NGC 4374 
NGC 4473 
NGC 4486 
NGC 4564 
NGC 4594 
NGC 4649 
NGC 4697 
NGC 4742 
NGC 5128 
NGC 5252 
NGC 5845 
NGC 6251 
Cygnus A 



Type 

S 
S 

s 
s 
s 
s 

E 
E 

s 

E 
S 
E 
E 
E 
E 
E 
S 
S 
E 
E 
E 
S 

s 

E 
E 
E 



Dist. 
[Mpc] 

0.8 
11.4 

7.5 

3.9 

9.7 
20.9 
11.2 
10.6 
11.6 
22.9 

7.2 
31.6 
26.2 
18.4 
15.7 
16.1 
15.0 

9.8 
16.8 
11.7 
15.5 

4.2 
94.4 
25.9 

101 
240 



M B 
[mag] 
-15.3 
-18.4 
-17.3 
-18.2 
-20.2 
-19.6 
-19.0 
-19.9 
-19.0 
-19.9 
-17.2 
-21.1 
-19.6 
-21.4 
-19.9 
-21.5 
-17.4 
-21.3 
-21.3 
-20.2 
-18.9 
-20.8 
-20.7 
-18.7 
-21.4 
-21.8 



M K 
[mag] 
-19.3 
-23.5 
-21.3 
-24.1 
-24.4 
-23.3 
-23.6 
-24.2 
-22.6 
-24.1 
-22.4 
-25.6 
-23.9 
-25.7 
-23.8 
-25.6 
-21.9 
-25.4 
-25.8 
-24.6 
-23.0 
-24.5 
-25.5 
-23.0 
-26.5 
-27.2 
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Figure 3. The long solid line shows the Mbh-L relation (equa- 
tion 1131 using our updated values for the galaxies in Marconi 
& Hunt (2003). The short lines emanating from the data points 
show the location of the galaxies as used by Marconi & Hunt to 
obtain the long dashed line. The location of IC 1459 as used by 
Marconi &; Hunt but excluded by us is shown by the cross. 



4-2.2 Removing questionable data points 

We repeated the above regression analysis removing four 
galaxies whose parameters are somewhat questionable. 
First, because the mass estimate for NGC 4742 has not 
yet appeared in a refereed paper, we hold off on its inclu- 
sion here. We also excluded NGC 1023 (r h = 0".08) and 
NGC 3377 (r h = 0".46) because their SMBH spheres of 
influence (rh, Merritt & Ferrarese 2001c) have apparently 
not been resolved according to Table II in Ferrarese & Ford 
(2005) in which rh/r res ratios of 0.89 and 0.74 respectively. 
Finally, as noted by Ferrarese & Ford (2005), NGC 4594 
has not yet had its SMBH mass acquired using a 3-integral 
model, and so it mass may therefore be in error. With our 
reduced sample of 22 objects, we obtain the relation 

log(M bh /M ) = -0.37(±0.04)[Mx+24]+8.29(±0.08), (14) 

consistent with the expression given in equation 1131 The 
total scatter about this relation is 0.33 dex, and the intrinsic 
scatter is O.SOt^os dex in logM bh . 

Comparison of equation ll4l (and ll3[) with our revised es- 
timate of McLure & Dunlop's /•f-band relation (equation [SJ 
reveals that we have resolved the disagreement noted in our 
Introduction. That is, our updated data sets and reanalysis 
of the McLure & Dunlop and the Marconi & Hunt stud- 
ies has yielded Mbh-L relations that agree with each other. 
Our preference is to use equation [14] because a) it was de- 
rived using both elliptical and disc galaxies and b) it has the 
smallest uncertainty on the slope and intercept. 



4-2.3 Is the Afbh-i relation non-linear 

Using the above 22 data points, we have explored whether 
the Mbh-i relation may be curved. The optimal log- 
quadratic relation, fitted in the same way as the Mbh-n data 
in Graham & Driver (2007a), is 



log(M bh /M ) = 



-0.37(±0.05) [M K +24] 
-0.01(±0.02)[Mk +24] 2 



8.32(±0.10) 
(15) 



The coefficient in front of the quadratic term is consistent 
with a value of zero, or in other words, it does not deviate 
from a value of zero at even the 1-sigma level. The 3a range 
of values on this term is only ±0.05. One can therefore con- 
clude that the Mbh-L relation, defined with the present data 
set, is not curved. 



4-2.4 Is the Mbh-L relation the same for elliptical 
galaxies and bulges 

For the 12 elliptical galaxies which comprise the sample of 
22 galaxies used in Section \A. 2. 21 the best-fitting relation is 

log(A/ bh /M ) = -0.33(±0.09)[Af A -+24]+8.33(±0.15).(16) 

For the 10 disc galaxies, which includes NGC 221 and 
NGC 4564, one has 

log(Af bh /M ) = -0.39(±0.08)[M A -+24]+8.33(±0.16).(17) 

Consequently, there appears to be no significant difference 
between the relations defined by the elliptical galaxies and 
the bulges of disc galaxies. 
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4-2.5 An R-band Mhh-L relation 

The above sample of 22 galaxies consists of 12 E galaxies, 
8 SO galaxies, and only one Sb and one Sbc galaxy. Using 
an R c — K colour of 2.6 and 2.5 for the elliptical and lentic- 
ular galaxies respectively, and 2.3 for the late-type galaxies 
(Buzzoni 2005) we obtained the following .R-band Mhh-L 
relation 

log(M b h/Af Q ) = -0.38(±0.04)[M K +21]+8.12(±0.08), (18) 

While this relation has overlapping error bars with the two 
independent R-band relations from the previous sections 
(equation[7]and[T2|, it is preferred for the reasons mentioned 
in Section 14^21 

4-2.6 A B-band Mhh-L relation 

Starting with the B-band magnitudes tabulated in Marconi 
& Hunt, (which have predominantly come from Tremainc 
et al. 2002 via Faber et al. 1997), we have modified these 
according to Section [4727T] Adjusting also the SMBH masses 
as in Section l4.2.1l we obtain the updated B-band expression 

log(M bh /M ) = -0.40(±0.05)[M s +19.5]+8.27(±0.08), (19) 

with a total scatter of 0.34 dex in logMbh- Using M B q = 
5.47 mag (Cox 2000), equation[l9]has a marginally shallower 
slope than the solution in Marconi & Hunt: log(Mbh/AfQ) = 
-0.48(±0.05)[M S + 19.53] + 8.18(±0.08). 

This new B-band slope of —0.40 ± 0.05 (obtained us- 
ing 22 galaxies) is comparable to the slope —0.42 ± 0.06 
(equation [5]) obtained in Section [3.11 using M32 plus the 16 
elliptical galaxies from McLure & Dunlop (2002). However, 
given no correction for dust attenuation was used, and given 
that the past B-band bulge/disc separation could probably 
be improved upon, this may perhaps be fortuitous. 

4-2.7 Scatter in the Mbh-Af sp heroid relation 

Marconi & Hunt (2003) additionally presented a relation 
between the mass of the black hole and the mass of the 
host spheroid. They used Af sp hcroid = 3R c o-^/G, in which 
Be and a e are the effective half-light radius and velocity 
dispersion of each spheroid, respectively. They reported (for 
their 27 Group 1 galaxies) an intrinsic scatter of 0.25 dex 
(and 0.49 dex for the full sample). Using their Group 1 data, 
and fitting equation [2] — which is designed to minimise the 
scatter in log Afbh — we measure the total scatter to be 
0.30 dex. A fuller, proper investigation of this subject would 
however require checking the R c values, which is beyond the 
scope of the current paper. 



5 DISCUSSION AND CONCLUSIONS 

Surprisingly, given the recent large body of work on SMBHs, 
there remains a strong need for high quality images for the 
sample of inactive galaxies with direct measurements of their 
supermassive black hole masses. In particular, near-infrared 
NICMOS images would enable an analysis of the core struc- 
ture of these galaxies. Graham & Driver (2007a) have pro- 
posed that the central stellar density may be the fundamen- 
tal parameter related to the mass of the black hole. This 



contrasts with current studies, including this one, which ex- 
plore the SMBH connection with global rather than nuclear 
properties of the host spheroid. Given the known trend be- 
tween central stellar density and host spheroid luminosity 
(e.g., Graham & Guzman 2003; Merritt 2006, his Fig.5), 
the popular relations may all be secondary in nature, i.e., 
subsequential Specifically, the central density (prior to core- 
depletiorjf] in massive spheroids — requiring the use of the 
core-Sersic model, Graham et al. 2003), or the central den- 
sity of less massive spheroids after modelling and exclud- 
ing the flux from their additional nuclear components (e.g. 
Graham & Guzman 2003) may be the key parameter con- 
nected to the SMBH mass. Due to the need for a bulge/disc 
decomposition for half of the current galaxy sample, such 
HST images should be mated with deep larger field-of-view 
ground-based images (e.g. Balcells et al. 2004) so as to ade- 
quately sample the domain of the disc. 

Aside from the study by Erwin et al. (2004) with 13 
galaxies, all of the optical Mhh-L relations to date have 
been constructed using the aperture growth curve magni- 
tudes obtained some 20 years ago and presented in Faber 
et al. (1989) or from the RC3 (de Vaucouleurs et al. 1991). 
A homogeneous set of high-resolution, deep, wide-field CCD 
images in the optical bands such as B and R would be highly 
useful for a) properly calibrating the Mhh-L relation, b) ac- 
quiring accurate bulge sizes and central densities (both pro- 
jected and deprojected) and c) subsequently calibrating the 
Afbh- (spheroidal mass) relation (e.g. Marconi & Hunt 2003; 
Haring & Rix 2004). Not only would this allow a proper 
test for the optimal fundamental relation, but it would pro- 
vide the community with improved relations for predicting 
SMBH masses in other galaxies. 

Nonetheless, there does now appear to be agreement 
between the various Mhh-L relations (see Table |3J. 

We plan to apply our B-band Mhh-L relation (equa- 
tion [T9j) to the Millennium Galaxy Catalogue (e.g., Driver 
ct al. 2006). This catalogue contains structural parameters 
from the B 1/ ^ n -bulge plus exponential-disc decomposition of 
10,095 nearby (z ~ 0.1) galaxies (Allen et al. 2006). While 
we have already been able to use the Afbh-n relation in Gra- 
ham & Driver (2007a) to predict the SMBH masses in these 
galaxies (Graham et al. 2007), the inconsistencies in the pre- 
viously published Mhh-L relations had prohibited their use. 
This will allow us to construct an updated SMBH mass func- 
tion for both early- and late-type galaxies, which can then 
be integrated to obtain the local SMBH mass density. 

5.1 Intermediate mass black holes 

Using the linewidth-luminosity-mass scaling relation by 
Greene & Ho (2005), Dong et al. (2007) report on the exis- 
tence of a 7 x 10 4 Mq intermediate mass black hole (IMBH) 
in the dwarf disk galaxy SDSS J160531. 84+174826.1 (see 
Figure 0}. Applying McLure & Dunlop's (2002) B-band 
Mhh-L relation to the central bulge/bar magnitude of this 
disk galaxy, they obtained a black hole mass one order of 

9 The apparent depletion of stars at the centres of giant galaxies 
may not have (only) arisen from the scouring action of coalescing 
SMBHs (see, for example, Boylan-Kolchin et al. 2004 and Nipoti 
ct al. 2006). 
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Table 4. New M^-L relations for predicting SMBH masses. The "origin" of the data is as 
follows. KG 2001 = Kormendy & Gebhardt (2001), not modified. MD 2002 = McLure & Dunlop 
(2002), modified (see Section |3T] and Table QJ. EGC 2004 = Erwin et al. (2004), modified (see 
Section |4~T1 and Table [2}. MH 2003 = Marconi & Hunt (2003), modified (see Section gj). and 
Tablc[3]l- The "sample" may consist of elliptical (E) galaxies or disc galaxies, denoted by 'S' for 
either SO or Sp. The total scatter in the log M bh direction is denoted A to t , while the intrinsic 
scatter in the logM D h direction is denoted 6i n t r i n sic- 
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Figure 4. SDSS J160531. 844-174826.1 (cross) has been added 
to the data points used to construct equation [7] shown here by 
the solid line. This dwarf galaxy, from Dong et al. (2007), has an 
intermediate mass black hole that has not been used in the fitting 
of the lines shown. The dashed line is equation I18l and the dotted 
line is equation [12] 



magnitude smaller than obtained with Marconi & Hunt's 

(2003) relation. Dong et al. give a (Ho = 70 km s _1 Mpc -1 ) 
i?-band magnitude of —13.9 mag. Adjusting this by 0.1 mag, 
to match our adopted value of -Ho = 73 km s _1 Mpc -1 , our 
updated i?-band expression from equation [7] predicts a black 
hole mass of 2.4 x 10 s Mq, while our updated relation in 
equation 1181 also gives a value of 2.4 x 10 5 Mq. This com- 
pares well with the value of 1.5 x 10 5 Mq obtained using the 
Mbh-n relation from Graham & Driver (2007a). 

We do not include the Seyfert 1 galaxy POX 52 as it 
appears to require a (yet to be performed) bulge/disc de- 
composition, with the disk accounting for the excess flux 
from 4-10 arcseconds seen in the Sersic fit of Barth et al. 

(2004) . Ideally, a bulge/disc decomposition is also required 
for the dwarf Seyfert 1, Sd galaxy NGC 4395. Inspection 
of figure 3 in Filippenko & Ho (2003) suggests that the ex- 
cess flux seen from 1 to ~2 arcseconds, peaking at ~0.6 
mag arcsec -2 above their fitted model, may be the bulge 
component of this late-typpe galaxy. However, due to the 



non-thermal point-source emission and nuclear-star cluster 
in this galaxy, it is difficult to know what may be the bulge 
component. While it appears that a Gaussian-like compo- 
nent should work well for the suggested bulge, we make no 
attempt to undertake this task here, but highlight the value 
of such a future investigation. The addition of more data 
points at the low-mass end of the M^-L relation should 
provide further valuable clues as to a) the range and relia- 
bility of this relation, and b) help constrain models for the 
co-evolution of black holes and galaxies. 

5.2 Implications for SMBH-galaxy coevolution 

Before discussing intrinsic physical relations, we need to per- 
form a symmetrical regression of Mbh and the magnitude. 
We have already seen that regressing Mbh on Mk gives a 
K-b&nA slope of —0.37 ± 0.04 (equation 1141) . Reversing the 
regression (see the text after equation [2} gives the relation 
log(M bh /M ) = -0.44(±0.05)[Af x + 24] + 8.30(±0.09), and 
so the average slope is ~ 0.40. This is consistent with the 
slope of — 0.454:0.05 from the symmetrical regression in Mar- 
coni & Hunt (2003). It is interesting to note that a slope 
of —0.40 implies a linear scaling between luminosity and 
black hole mass, such that M bh cc L 1()0 . That is, the SMBH 
mass to spheroid luminosity ratio is constant, as predicted 
by Chien (2007). 

Our updated data and new relation supports the "dry 
merger" scenario in which the SMBH mass doubles as the 
luminosity doubles. While this paints a consistent picture 
at the high-mass end, where "core galaxies" exist and dry 
mergers are considered the order of the day, this is unlikely 
to hold at lower masses. 

SMBH accretion of ISM gas derived from stellar winds 
is believed to fuel some AGN activity today (Fabbiano et 
al. 2004, Pellegrini 2005; Soria et al. 2006). Indeed, stellar 
mass loss in elliptical galaxies produces > 10 2 times more 
mass than that found in the central SMBH (Ciotti & Os- 
triker 2007). As these Authors note, when fuelling SMBHs 
with the recycled gas from the stellar population, almost 
by definition, the amount of fuel is proportional to the host 
galaxy's stellar mass. Ciotti & Ostriker also show how radia- 
tive heating from such AGN feedback is responsible for the 
self-regulated coevolution of galaxy and SMBH. Our results 
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provide a valuable new constraint on these models, setting 
the proportionality constant to one. 

Moreover, the level of total scatter in the Mbh-L relation 
(0.31 to 0.34 dex) makes it competitive with both the M^-n 
(Graham et al. 2001; Graham & Driver 2007a: 0.31 dex) and 
the Mbh-o" relation (Ferrarese & Merritt 2000; Gebhardt et 
al. 2000). While modelling the 31 galaxies from Tremaine 
et al. (2002) gives an intrinsic scatter of 0.27 dex, the total 
scatter about the relation presented there is 0.34 dex. 

This opens the question as to whether the stellar (bary- 
onic) or total (stars plus dark matter) mass is the driver 
of the SMBH-galaxy connection. Which of these, or some 
other quantity may be the fundamental parameter connect- 
ing galaxies and their black holes is not yet clear. 
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APPENDIX A: POTENTIAL SAMPLE BIAS 

Bernardi et al. (2007) have suggested that the local sample 
of inactive galaxies with direct SMBH mass measurements 
may be biased with regard to the greater population. They 
argue that, relative to the total population, there is a bias 
such that these galaxies have overly large velocity disper- 
sions for their luminosities. If correct, the result is such that 
the Mbh~i relation will over-predict the SMBH masses (in 
other galaxies) relative to the M^-o relation. We explore 
this issue with our updated data set. 

Using r - R = 0.24 (Fukugita 1995), the SDSS L-o 
relation from Tundo et al. (2007, their equation 4) is such 
that 

logo- = 0.27 - 0.092M fl . (Al) 

We compare this with our updated i?-band data from Erwin 
et al. (2004) and our updated if-band data from Marconi 
& Hunt transformed into the 7?-band using the tables in 
Buzzoni (2005), as done in Section 14.2.51 The results are 
shown in Figure IA1I 

Applying equation [2] to our updated Marconi & Hunt 
data set, in which we have minimised the scatter in the log o 
direction, we obtain 

logo- = 2.268 - 0.082(M fl + 21), (A2) 

which is shown by the solid line in Figure IA1I Compared to 
the dashed line (from the SDSS data set), the local sample 
of inactive galaxies therefore appears to have larger velocity 
dispersions for a given magnitude. 

While the local sample of galaxies with direct SMBH 
mass measurements may be biased, it is also possible that 
the magnitudes used to construct the SDSS L-o relation may 
have been under-estimated due to dust, or over-estimated 
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Figure Al. The solid line (equation \k2\ shows the regression of 
log <r on Mr using our updated Marconi & Hunt data set, denoted 
by the open squares. Our updated Erwin et al. data points are 
shown by the filled circles but are not used in the linear regres- 
sion shown here. The dashed line (equation |AH is the equivalent 
regression using SDSS data (Tundo et al. 2007, their equation 4). 
A slight offset is evident. 

due to the bulge-disc separation. That is, they may be the 
biased data set. Unlike the Erwin et al. (2004) and Mar- 
coni & Hunt (2003) bulge data, no Sersic i? 1/n -bulge + 
exponential-disc fit was performed in acquiring the SDSS 
bulge magnitudes. The use of R 1 ^ models is known to over- 
estimate the bulge flux for bulges with stellar distributions 
having n < 4 (e.g. Graham & Driver 2007b, their Table 3), 
which may be the bulk of the SO galaxy population (e.g., 
Balcells et al. 2003) . A second concern is that faint elliptical 
galaxies have a different slope in the L-o diagram than lu- 
minous elliptical galaxies. The successive inclusion of fainter 
galaxies should therefore lead to a progressive change in the 
slope of the L-o relation. 

Ignoring the above issues for now, the average veloc- 
ity dispersion excess is given by differencing the above two 
equations, to give A logo- = 0.30 + 0.011M.R. Therefore, if 
one wished to correct for this (possible) sample bias, the 
local galaxy magnitudes Mr (i.e., those with direct SMBH 
mass measurements) could be adjusted by AM, such that 

AM = 0.109[M« + 27.60]. (A3) 

This would bring the two lines in Figure |AH into agreement. 
The same trick was performed in Tundo et al. (2007, at 
the end of their Section 3), although using different bulge 
magnitudes to our updated values. In a future paper we 
intend to better address such a potential offset in the L- 
o diagram using the dust-corrected ii 1 ^ n -bulge magnitudes 
from the Millennium Galaxy Catalogue's 10 095 galaxies 
(Allen et al. 2006). 

At the low luminosity end the situation is different to 
presented above. This is because the L-o relation is not lin- 
ear, having a well-recognised slope of ~4 at the bright end 
and a less well known, but long established, slope of ~2 
at the faint end (Tonry 1981; Davies et al. 1983; Held et al. 
1992; De Rijcke et al. 2005). Matkovic & Guzman (2005) ar- 
gue that the transitional magnitude occurs near Mr ~ —22 
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Figure A2. The solid and dashed lines have the same meaning 
as in Figure IA1I as do the squares and large filled circles. The 
small dots are data taken from Matkovic & Guzman (2005, their 
Figure 4, adjusted from H =70 to 73 km s _1 Mpc -1 ) and the 
thick solid line is their /Zo-adjusted (logo - | Mr) regression for 
faint galaxies. 

mag, coinciding with the onset of dry merging and the break 
seen in the luminosity - central surface brightness diagram 
(Graham & Guzman 2003, their figure 9c). What is also ap- 
parent in Figure 4 from Matkovic & Guzman is the subtle 
nature of this transition, and the need for a long baseline in 
magnitude for it to be recognised. 

Fitting a single power-law to the L-a data should yield 
a slope that is dependent on one's luminosity range and 
thus sample selection. This has not received much attention 
in the literature, most likely because of the subtle nature 
of the transition due to the scatter about the relation. For 
this reason, the change in slope only really becomes obvious 
in samples containing magnitudes fainter than Mr ~ — 20 
mag. 

Performing a regression of log a on Mr, Matkovic & 
Guzman (2005) report a faint-end relation of 

Mr = (-5.585 ± 0.210) log a - 8.755(±0.444), (A4) 

or simply 

log cr = -0.179Mr - 1.508, (A5) 

which has been adjusted here to Ho=7S km s _1 Mpc - . 
Their data and this relation can be seen in Figure lA"2l along 
with our sample of galaxies with direct SMBH mass mea- 
surements. One can see that the low-luminosity extrapo- 
lation of the L-sigma relation obtained using the updated 
Marconi & Hunt data set (equation IA2|) does not follow 
the trend defined by the fainter galaxy population (equa- 
tion \KE§. The problem is such that for a given velocity dis- 
persion, the magnitudes predicted from the (SMBH sample)- 
derived L-a relation are too faint, compared to the general 
population. The magnitudes predicted using equation IA2I 
need to brightened by 

AM = 0.542[Mr + 21.78]. (A6) 

If the sample of galaxies with direct SMBH masses have 
normal masses and velocity dispersions, but biased lumi- 




Figure A3. The data points have the same meaning as in Fig- 
ure \3\ although here we only show those 22 galaxies from Sec- 
tion [47372] The short lines emanating from each data point show 
their location prior to the magnitude- adjustment performed in 
Section fS] The location of the data points show their magnitude- 
adjusted location. The dashed line shows the regression prior to 
this adjustment ( equation 1186 while the solid line shows the new 
regression ( equation IA7j l, 

nosities, as argued by Bernardi et al. (2007), then one can 
apply the above magnitude corrections to this sample. At 
Mr = —20.28 mag we switch from using corrective equa- 
tion's] to corrective equation I A6I With our modified set of 
galaxy magnitudes, corrected to represent the general pop- 
ulation, one can repeat the regression analysis to obtain a 
new Mbh-L relation. Doing so, one has 

log(Af b h/M ) = -0.51(±0.06)[Mr+22]+8.16(±0.09), (A7) 

with a total scatter of 0.36 dex, and an intrinsic scatter of 

u - o± -0.06- 

Compared to equation 1181 which was obtained prior 
to this magnitude adjustment, for magnitudes fainter than 
Mr ~ —24.5 mag the new relation predicts smaller SMBH 
masses (see Figure IA3[) . Extrapolation of the relations to 
brighter magnitudes results in greater SMBH masses. This 
is a consequence of the L-a relations from Figure IA1I that 
we have used to derive the magnitude adjustment in equa- 
tion [M}. 

An alternative scenario is that the M^-L relation may 
have two slopes, described by a broken power-law with the 
transition denoting the onset of dry merging. Indeed, the 
above prescription should generate such a relation. If one 
accepts that Mbh oc a 4 , then at the luminous end, where 
L cc a 4 , one naturally obtains Afbh oc L . For magnitudes 
fainter than Mr ~ —22 mag, one has L oc a 2 , and if Mbh oc 
a 4 over this domain, then one should expect to find Aibh oc 
L 5 . However, we feel that to properly address this scenario 
will require more data than is available at present. 



